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Abstract 

Positronium is an ideal system for the research of the bound state 
QED. New precise measurement of orthopositronium decay rate has 
been performed with an accuracy of 150 ppm. This result is consistent 
with the last three results and also the 2nd order correction. The 
result combined with the last three is 7.0401ib0.0007/xsec~^ (100 ppm), 
which is consistent with the 2nd order correction and differs from the 
1st order calculation by 2. 60". It is the first test to validate the 2nd 
order correction. 



1 Introduction 



Positronium (Ps), the bound state of an electron and a positron, is a purely 
leptonic system and the triplet (l^S*!) state of Ps, orthopositronium(o-Ps), 
decays slowly into three photons. Precise measurement of this decay rate 
gives us direct information about quantum electrodynamics(QED) in bound 
state. 

Three precise measurements [H El [3] of the o-PS decay rate were per- 
formed, in which reported decay rate values were much larger, i.e., 5.2 - 9.1 
experimental standard deviations, than a QED prediction [H [5] (7.0382 /isec~^) 
corrected up to 0{a) and 0(a^ log(l/a)). This discrepancy has been re- 
ferred as 'orthopositronium lifetime puzzle', and was long-standing problem. 
In an effort to clarify the discrepancy, various exotic decay modes have been 
searched for, without any evidence so farP [3, El El flO] . 

It was found [HI |12] that thermalization process of the produced o- 
Ps is much slower than its lifetime, and the unthermalized o-Ps becomes 
the serious systematic error in all the previous measurements. With new 
measurement [11]. we obtained decay rate of 7.0398(29) /isec~^, which was 
consistent with the first order QED calculation, and quite differ from the 
previous results. In this method, the thermalization process was directly 
measured and taken into account. This result has been confirmed by two dif- 
ferent and more accurate (200 ppm) measurements [T^ [T3] and the 'lifetime 
puzzle' has been solved p!5]. Since the same results are obtained with differ- 
ent two methods([lIl [13] [E]), these methods have the potential to measure 
decay rate more precisely. 

In our previous two measurements [TT| [13] . the Nal(Tl) and CsI(Tl) scin- 
tillators were used, whose decay constants were long and rise times were also 
not fast. It was the reason that the timing performance of the previous mea- 
surements were limited. Furthermore, a low activity source of the positron 
only can be used in order to suppress pileup, since they have long decay con- 
stants. Many types of newly developed scintillators have been tested, and 
the YAP scintillator is used for this measurement. We can obtain good time 
performance and we can use positron source with an activity of 11 kBq. In 
order to reduce the trigger of the useless events, anti-trigger system is intro- 
duced. If the positron annihilate into 27 on the vacuum vessel, trigger signal 
is vetoed. 

Non-relativistic-QED approximation (NRQED) has been developed re- 
cently, and it is useful to calculate the higher order correction of the bound 
state. This approach can be used for not only QED but also QCD, in which 
the higher order corrections are crucial due to the strong coupling constant. 
The second order correction (9(a^)|l6], whose contribution is about 243 ppm. 
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has been performed in 2000 and 2002, and the more accurate measurement 
is necessary to examine this higher order correction. The logarithmic cor- 
rections 0(q;^ log(l/Q;)) and 0(a^ log^(l/a)) are also calculated and these 
contributions are found to be small, 3-4 ppm[T7l[IH]- In this paper, we focus 
on the validation of the 0{a^) correction. 

2 Experimental setup 

Figure 1 shows a diagram of the experimental setup [I9]; a ^^Ge — Ga positron 
source with the strength of 11 kBq, being sandwiched between two sheets 
of plastic scintillators(NE102 thickness=200/im) and held by a cone made 
of aluminized mylar. The cone was situated at the center of a cylindrical 
vacuum container made of 1 mm-thick plastic scintillators and glass, being 
filled with Si02 aerogel (RUN-I) or powder (RUN-II), and evacuated down 
to 1 X 10^^ Torr. Density of the aerogel and powder is 0.03 g/cm^ both, 
and the surface of the primary grain are replaced into hydrophobic in order 
to remove the electric dipole of the 0H-. The sizes of the primary grain 
are 10 and 16 nm for the aerogel and powder, respectively. Thus the mean 
free path of the collision between positronium and Si02 are different, and it 
makes difference in the thermalization process and the pick-off probability. 

The emitted positrons pass through thin scintillator, transmitting a light 
pulse to 1/2-inch photomultiplier(Trigger PMT:Hamamatsu H3165P-10VT) 
and forming Ps when stopped in the silica powder or aerogel. The signal of 
the Trigger PMT is used for start signal(t=0) of all TDCs. Since the kinetic 
energy of the positron is high (ii^end= l-9MeV: ^^Ge — Ga) about 40% of the 
positron pass through the silica and annihilate at the vacuum container. If 
the positron deposit energy(>100 keV) in the 1 mm-thick plastic scintillators, 
inside wall of the vacuum container, the produced trigger signal is vetoed to 
remove these annihilate events in the wall. A light pulse from the 1 mm-thick 
plastic scintillators is guided into two 1-inch photomultipliers(Anti-trigger 
PMT:Hamamatsu R1924P), and gives the anti-trigger signal. 

Three high-purity coaxial germanium detectors (Ortec GEM38195) pre- 
cisely measured the thermalization process. One output from the pream- 
plifier of each germanium detector is fed through a timing-filter amplifier 
(Ortec 474), whose integration and differential times are set to 20 and 500 
nsec, respectively. Outputs of the timing-filter amplifier are put into three 
leading-edge discriminators, in which different threshold values are set at 30, 
50 and 150 keV, and these three timings are recorded individually by TDCs. 
The signal shape of the coaxial germanium detectors depends on the posi- 
tion at which energy is deposited, since the electric field for charge collection 
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Figure 1: Schematic diagram of apparatus. 



is complicated. Time walk correction can be applied and correct timing is 
obtained using these three timings. The good time resolution of (4 nsec) is 
obtained for 480-511 keV. 

Another output from the preamplifier of each germanium detector is fed 
into a spectroscopy amplifier (Ortec 673), and the signal is integrated with the 
time constant of 6/isec. The output is recorded with peak-hold ADC(Hoshin 
coil). Energy resolutions were measured using several line 7 sources, with 
typical resultant values of 0.5 keV at 514 keV. A 2 mm-thick lead sheet is 
placed in front of each germanium detector to eliminate the contribution of 
the two low-energy 7's from the three-photon decay of o-Ps. 
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Four YAP(YA103: Ce doped: size is 50x50x33 mm) scintillators si- 
multaneously measured the time and energy information from each decay. 

2- inch photomultiplier(Hamamatsu R329-02) is connected to each crystal. 
The charge ADCs are used to measure the energy information. A single 
leading-edge discriminator (threshold 70 keV) is used for the timing, and 
the time walk correction is applied as a function of the energy information. 
Since the decay constant of the YAP is fast (30 nsec) and the slow component 
of the scintillation is negligibly small, the event pileup probability is small 
and the time walk correction can be performed perfectly as a function of 
the deposited energy. The good time resolution of 1.2 nsec is obtained for 

>150 keV. YAP is an ideal scintillator for the timing measurements; high 
density (5.27g/cm^), high light yield (40% of Nal(Tl)), low contamination 
of the radio activities, the fast rise time(<10 nsec), and the fast decay con- 
stant (30 nsec). Detail studies of the various scintillators are summarized in 

m- 

A direct clock (2 GHz) count type TDC is used for both the germa- 
nium detectors and the YAP scintillators. The time resolution and the full 
range are 0.5 nsec and 32 /isec(16bit), respectively. The used clock is cali- 
brated with an accuracy within 1 ppm and stable (< 2.5 ppm within —10 - 
50°C). The integrated non-linearity is negligibly small and the differential 
non-linearity is 0.1% for 4bin. 

Data are taken with the CAMAC system for 2.3 months(RUN-I) and 3.1 
months(RUN-II), and the trigger rates before the anti-trigger are 10.8 and 
6.8 kHz, respectively. Total observed events of j3^ emission are 1.4 x 10^'^ and 
1.6 X 10^°. The energy calibration and time calibration has been performed 
for every 6 hours. Area temperature was maintained within ±0.5°C to ensure 
the stability of energy and time calibration. 

As some fraction of o-Ps inevitably results in 'pick-off' annihilation due 
to collisions with atomic electrons of the target material, the observed o-Ps 
decay rate Xobs is a sum of the intrinsic o-Ps decay rate As^ and the pick-off 
annihilation rate into 27's, Apjcfc, i-e., 

Ao6s(^) = + Xpickit). (1) 

Since the velocity of the Ps slows down slowly due to the elastic collision 
to the material, pick-off probability depends on the time after produced as 
shown in Figures. 2. The energy distribution of photons emitted from the 

3- body decay is continuous below the steep edge at 511 keV, whereas the 
pick-off annihilation is 2-body which produces a 511 keV monochromatic 
peak. Energy and timing information are simultaneously measured with 
high-energy resolution germanium detectors such that Xpickit) / X^-y can be 
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determined from the energy spectrum of the emitted photon. Once a precise 
thermahzation function is obtained, Xpick{t) will contain all information about 
the process. The population of o-Ps at time t, N(t) can be expressed as 

N{t) = iV^ exp (^-A3, /J (l + dt''^ . (2) 

Providing the ratio is determined as a function of time, the intrinsic de- 
cay rate of o-Ps, As^, can be directly obtained by fitting the observed time 
spectrum. 



3 Analysis 

The ratio Xpickif) / ^z-y is determined using the energy spectrum measured by 
the germanium detectors. The energy spectrum of the o-Ps decay sample, 
referred to as the o-Ps spectrum, is obtained by subtracting accidental con- 
tributions from the measured spectrum. The 37-decay continuum spectrum 
is estimated using Monte Carlo simulation in which the geometry and vari- 
ous material distributions are reproduced in detail [19j. For every simulated 
event, three photons are generated. Successive photoelectric, Compton, or 
Rayleigh scattering interactions of every photon are then followed through 
the materials until all photon energy is either deposited or escapes from the 
detectors. The response function of the detectors is determined based on the 
measured spectrum of monochromatic 7-rays emitted from ^^^Eu, ^^Sr, and 
^^'^Cs. These material and detector effects are taken into account precisely, 
and the Sj-spectrum is obtained. These are continuous distributions, and is 
normalized to the observed o-Ps spectrum with the ratio of event numbers 
within the region (480-505 keV). 

Figures 2(a) (b) show the o-Ps spectra and 8-7 spectra observed in RUN- 
I(Aerogel) and Il(powder), respectively. The pick-off 2-7 spectra, which are 
obtained from the o-Ps spectra after subtracting the 37-spectra, are also 
superimposed in the same figures. The \pickif) / X^-y can be calculated di- 
rectly from the ratio of the pick-off 2-7 and 8-7 spectrum, the ratio does 
not depend on the absolute values of the detection efficiencies. The relative 
value of the detection efficiency esu fcey/e480-505 keV is estimated with the 
monochromatic 7-rays emitted from ^^^Eu, ^^Sr, and ^^^Cs. The calculations 
of the \pick{t) / ^s-y are performed for the various time windows, and its time 
dependence are shown in Figure 2(c)(RUN-I) and (d)(RUN-II). Horizontal 
axis of the figures is the time between the positrons emission and decay, this 
slope shows the thermahzation process of the positronium. It takes much 
time than the lifetime to well thermalized, as we have already shown in the 
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previous measurements pT| [T3j . Since the collision rate in RUN-I is higher 
than that in RUN-II, slopes of the thermalization process are different as 
shown in these figures. 




Figure 2: (a) Energy spectra of o-Ps decay 7's obtained by germanium de- 
tectors in RUN-I(Aerogel). The solid line represents data points in a time 
window of 60-700 nsec, and the dotted line shows the 37-decay spectrum 
calculated by the Monte Carlo simulation. Pick-off spectrum obtained after 
subtracting the 87 contribution from the o-Ps spectrum is superimposed, 
(c) The ratio Xpickit) / ^z-y are plotted as a function of time (RUN-I). Only 
the statistical errors are shown and the solid lines represent best fit results 
obtained. (b),(d) the same figures for RUN-II (Powder) 

Figures 3(a) and (b) show the time spectra measured with the YAP scin- 
tillators in RUN-I and -II. A sharp peak of the prompt annihilation is followed 
by the exponential decay curve of o-Ps and then the constant accidental spec- 
trum. Time resolution of 1.2 nsec is obtained and the o-Ps curve is widely 
observed over 1.0 /isec. We fit resultant time spectrum using the least square 
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method with the following function; 



Nobs{t) = e^{-Rstopt) 

where Spick and are the detection efficiencies for pick-off annihilation and 
87 decays, and Rstop is an experimental random counting rate representing 
the fact that time interval measurement always accept the first 7 as a stop 
signal. Xpickit) / ^s-y is about 1% due to the low-density of the Si02 powder or 
aerogel. This means that the ratio of error propagation to the decay rate is 
suppressed by a factor of 100. It is a merit of our method, and an accuracy 
of 100 ppm can be achieved if Xpick{t) / ^S'y is measured with an accuracy of 
1%, which is the reafistic. 



l + ^^^^l N{t) + C 



^37 



A 



37 



(3) 
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Figure 3: (a) Timing spectrum measured by the YAP scintillators within an 
energy window above 150 keV (RUN-I). The figures in the right corner show 

the magnified view of the prompt peak, (b) The same timing spectrum in 
RUN-Il. (c) (d) The fitted results are plotted as a function of the fitting start 
time for RUN-1 and -11. 

The fitting procedure has been performed for the various fitting region. 
Figures 3(c) and (d) show the fitted results as a function of the fitting start 
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Source 


RUN-I(ppm) 


RUN-II(ppm) 


TDC noil linearity 


< ± 15 


< ± 15 


Pile iin 


< + 10 


< + 10 


JT IL^lvUii collllidllUii 






ill '^'"i/ en rTf" TSl (^"i'lOTl 
I ± I ^ / oU. U LI Llwli 


\ _1_ Oc/ 


^ + Q1 


(2) Ge efficiency 


< ± 33 


< ± 28 


(3) YAP efficiency 


< ± 64 


< ± 19 


Physics source 






(1) Zeeman effect 


< - 5 


< - 5 


(2) Stark effect 


< + 3 


< + 4 


(3) 37 annihilation 


< - 91 


< - 33 


Total 


-147 and + 115 


-104 and +98 



Table 1: Summary of the systematic errors 

time, end point is fixed at t=3600 nsec. If the ratio Xpickit) / ^3-/ would be 
mismeasured, the ffited results becomes unstable against the fitting region. 
The fitted results are stable against the fitting regions as shown in these 
figures and the values with the fitting start time of 60 nsec are used for both 
run. The obtained decay rates[T9j are As^ = 7.03876±0.0009(stat.) /isec~^ for 
RUN-I and 7.04136 ± 0.0009(stat.) fisec'^ for RUN-II, which are consistent 
with each other. 

Estimates of various systematic errors are summarized in Table 1 . 



The predominant contribution to total systematic error is produced by 
uncertain normalization of the 3-7 spectrum. That is, the number of pick-off 
events are determined by subtracting the normalized 3-7 spectrum of Monte 
Carlo simulation from the o-Ps spectrum, where changing the normalization 
factor affects the Xpickif) / A37 values and eventually propagates to the final re- 
sult. Since the sharp fall-off of the 37-spectrum at 511 keV is solely produced 
by the good germanium energy resolution of a = 0.5 keV, this subtraction 
only affects the lower side of the pick-off spectrum such that improper sub- 
traction results in asymmetry of the pick-off spectrum shape. Comparison 
of the asymmetries of the pick-off peak shape and the prompt peak annihila- 
tion spectrum is a good parameter for estimating this systematic error. The 
errors are about 90 ppm for both measurements. 

The relative value of the detection efficiency esu fcey /e480-505 keV is used 
to estimate Xpickif) I X^-y- The relative efficiency is directly measured using 
the 7-rays and there are two other uncertainties in the relative efficiency. (1) 
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If the density of silica power(RUN-II) changes by 10% during the run time of 
3 months, the distribution of the decay points of orthopositronium changes 
and detection efficiency change by 0.12%. (2) Some part of the germanium 
signal has very slow rise time(SRT), and these events are rejected using the 
timing information of three different thresholds. This SRT rejection efficiency 
are estimated using the real data. The uncertainties of the estimation are 
0.08%(RUN-I) and 0.09%(RUN-II). These two uncertainties are propagated 
to the fitted valeus of As^, about 30 ppm for both runs. 

The fitting function includes the relative efficiency of the YAP scintillator 
as shown in the equation(3). This relative efficiencies are also estimated with 
real data using the 7-rays emitted from ^^Sr(514 keV). An uncertainties of the 
relative efficiency are 2%(RUN-I) and 1.3%(RUN-II), and they contribute to 
the errors on the fitted valeus of X^^ by 64 ppm(RUN-I) and 19 ppm(RUN-II). 

The Stark shift stretches the lifetime of Ps atoms. A calculation shows 
that the shift is proportional to a square of the effective electric field E such 
that AAa^/Ag^ = 248 • {E/EqY, where Eq = mle^h^ ^ 5.14 x 10^ V/cm. 
Silanol functional groups on the surface of the grain behave as an electrical 
dipole moment creating an effective field around the grains. Average den- 
sities of Silanol are measured to be 0.44nm~^ and the effective field can be 
analytically calculated such that the contribution to the o-Ps decay rate is de- 
termined to be -5 ppm. These estimations were confirmed using results from 
precise hyper-fine-structure (HPS) interval measurements of ground state Ps 
in silica powder [20], where the interval is proportional to the size of Stark 
effect. Considering the difference in powder densities used, the HPS results 
are consistent with our estimation. 

Error contribution due to the Zeeman effect is estimated using the mea- 
sured absolute magnetic field around the positronium assembly (-5 ppm). 
Since the 3-7 pick-off process can only occur at a certain ratio, and the cal- 
culated relative frequency. On the other hand, the pick-off events by the 
spin-flip effect do not have the 37 decay process. If the contribution of the 
spin-flip is neglected, the 37 pick-off events increase by 0.3 % and, being 
-91 ppm for RUN-I and -33 ppm for RUN-II. 

The above discussed systematic errors are regarded as independent con- 
tributions such that the total systematic error can be calculated as their 
quadratic sum, resulting in -147 ppm, -|-115 ppm for RUN-I and -104 ppm, 
-f98 ppm for RUN-II. The combined result with systematic error is As^ = 
7.0401 ± OmOQ{stat.)toMolisys.) ^sec-\ and total error is 150 ppm. Cor- 
relation of the systematic errors are taken into account. 
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4 Results and discussion 



Figure 4 shows the summary of the measured decay rate after 1995, three 
results are obtained with our method and one result is obtained using thin 
polymer, with which the produced positronium is very slow and the effect of 
the thermalized positronium is suppressed. These four results are consistent 
with each other. And the combined value is As^ = 7.0401 ± 0.0007 /isec~^ 
and shown in the red arrow in the figure. The correlations of the systematic 
errors are carefully taken into account. The combined result is consistent 
with the all four results and accuracy is 100 ppm. This result is consistent 
with the 0(a;^) correction ^16j. and differs from the only up to 0{a) by 2.6a. 



without 0(a' 
correction 



QED prediction 
witli O(a^) correction 



Phys.Lett. B357,475 (1995) 



Phys.Rev.Lett. 90, 203402 (2C 02) 



Phys.Lett. 6572,117(2003) 



This Measurement 



-► combined result 



7.037 7.038 7.039 7.04 7.041 7.042 7.043 

Decay rate (ii sec-^) 



Figure 4: Decay rates measured in the last four experiments are listed chron- 
ically. The last arrow shows the combined result of the last four measure- 
ments. The red and blue lines shows the QED prediction calculated up to 
0(a;^) and 0(a), respectively 
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5 Conclusion 



The decay rate of o-Ps was measured using a direct 27 correction method 
in which the thermahzation effect of o-Ps is accounted for and integrated 
into the time spectrum fitting procedure. The new result is As^ = 7.0401 ± 
Om06{stat.)tomolisys.) fisec'^, and total error is 150 ppm[in]. It is the 
most accurate result and agrees well with the last three results[TTl [T31 E]. 
The combined result of these four measurements is 7.0401±0.0007/isec~^, 
which is consistent well with the NRQED prediction corrected up to O(a^) 
term[T6], and differs from the result up to 0{a) by 2.6a. It is the first result 
to validate the O(a^) correction. 

We thank Prof. M.Kobayashi(KEK) for very useful discussion about the 
scintillators. Sincere gratitude is extend to Dr.M.Ikeno(KEK) and Dr. T.Namba(U. Tokyo) 
for a development of 2 GHz TDC and the great support to this measurement. 
We are pleased to acknowledge the Japan Society for the Promotion of Sci- 
ence for the financial support. 
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